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ABSTRACT 

An improved model for the Earth’s gravity field, TEG-1, has been determined using data sets 
from fourteen satellites, spanning the inclination ranges from 15° to 115°, and global surface gravity 
anomaly data. The satellite measurements include laser ranging data, doppler range-rate data, and 
satellite-to-ocean radar altimeter data measurements, which include the direct height measurement 
and the differenced measurements at ground track crossings (crossover measurements). Also 
determined was another gravity field model, TEG-1S, which included all the data sets in TEG-1 with 
the exception of direct altimeter data. The effort has included an intense scrutiny of the gravity field 
solution methodology. The estimated parameters included geopotential coefficients complete to 
degree and order 50 with selected higher order coefficients, ocean and solid Earth tide parameters, 
doppler tracking station coordinates and the quasi-stationary sea surface topography. Extensive error 
analysis and calibration of the formal covariance matrix indicate that the gravity field model is a 
significant improvement over previous models and can be used for general applications in geodesy. 

1. INTRODUCTION 

Significant progress has been achieved during the last decade in the determination of the spherical 
harmonic coefficients of the Earth’s external gravitational potential. A substantial portion of this 
progress can be directly attributed to the advent of Earth-orbiting artificial satellites and to the ability 
to observe their motion from either ground-based or satellite-originated tracking data. While the 
satellite data primarily resolve the long and intermediate wavelengths (> 1500 km), global surface 
gravity measurements and the altimeter data are capable of recovering the shorter wavelength 
components of the Earth’s gravity field. Recent trends in gravity model improvement have been 
driven, in part, by requirements for more accurate satellite orbits to achieve the objectives of the 
Crustal Dynamics Project and the recently approved NASA/CNES Topex/Poseidon mission. A joint 
effort to develop an improved model for the Earth’s gravity field has been undertaken to develop a 
gravity model to meet the orbit accuracy requirement of the Topex/Poseidon mission. The gravity 
field solution will represent the first complete reiteration of the historical tracking data used to define 
the NASA Goddard Space Flight Center (GSFC) Earth Model series. The GSFC GEM-T1 [Marsh et 
al, 1988] and the University of Texas (UT) TEG-1 fields, described in this paper, are preliminary 
versions of the Topex gravity field solution. 

2. THEORY AND METHOD 

The gravitational potential, U , due to the Earth’s nonspherical mass distribution can be expressed 
as follows 
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where GM is the product of the gravitational constant and the total mass of the Earth and the 
atmosphere; R e is the mean equatorial radius of^he Earth; P ^ are the normalized Legendre 
associated function of degree / and order m; C lm ,S lm are the the normalized spherical harmonic 
coefficients whose values are functions of the mass distribution within the Earth and the atmosphere; 
AC /m , AS im are the time-varying components of C lm and S !m caused by tides; also are functions of the 
tidal coefficients, C,± and S,±; and r,<D,X are the Earth-fixed spherical coordinate system; r is the 
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men^ian StanCe ’ * “ ** 8e ° centric latitude and X is longitude measured from the Greenwich 


e estimation of C/ m , S lm , C/ m ,S im and other orbit and geophysical parameters can be 

H " r ne a m0d,fie ? procedure. This esttmation procedure whi* 

p ovides adjustments to satellite orbit-dependent parameters and other geophysical and geodetic 

rel!T! ’ y y ,€y |19751 and mod ' fied “ «• simultaneous estimation Tto 

weights for the individual satellite information arrays [Yuan etal., 1988]: 

X = (H t R~ l H)~ l H t R~ x y ; /?,. = x, k . (y. _ H . iy \ (y. _ H . £) . J 

where jc is the state parameter; R is the weighting matrix; / is the identity matrix; //, is the partial 
denvattve wtth respect to a for the ith data set; and k, is the number of observations for ith darnel. 

rech!llr! 3 e 2 r d , 0 <) 71 ? iV ^ ab0Ve Ca ” b ' SOlVed “ en “ ively “ sin * »"h°8onal tnntsformation 
T 7716 estimation process has been implemented in the University of 

Texas Orbit Processor (UTOPIA) software system [Sckun and Tapley 19801. Hie rln Zl, 

Linear ^stL^olv^fLLi™ T " 0nSa ‘ ellite info ™ aIi <>" «lualions was installed in the Large 
Linear System Solver (LLISS). Vectorized versions of UTOPIA and LLISS are operational on the 

mversi y o exas System Center for High Performance Computing Cray X-MP/24 supercomputer 
e erence orbits for each of these satellites were computed using UTOPIA with the best a priori 
gravity model and gravity field information equations were generated for eth data set TOe 
combination solution was performed using LLISS. 6 

3. DATA AND MODELS 

characrerhnics ^nd'dara Z ” SeleC,ei1 ^ ** Snlvity m0del «»«*>". Their orbital 

Ja^Sm 5- to “ Tab " '■ 71,e of these satellite orbt.s 

>, . _ 0 Peo e to 115 for Geos-3. The solution includes data at 90° for Oscar- 14 and 

surfl Da ^ types 1 i ndude laser ran8e * one ' wa y range-rate, altimeter, altimeter cro^ov ^d 
^trfacegravtty data. Detatled descriptions of the gradational and nongrav adonal force3els to 

Eatth orrentahon and o me model, laser, doppler, direct altimeter and surface gradre toa^remem 
models are summarized in Tapley et al. [1987], gravity measurement 

4. SOLUTION 

The list of parameters which are simultaneously estimated with a relative weighting factor for 
each data set include: (1) geopotential complete to degree and onlerSO. plus selected coefficients (2) 

ddei g ™ ”? 3 P ," 8 Ttr°? <m = °' ' = 2 ’ 3): Saa > Sa ' 

K2 and T2 cl - 2v' rdl'n ' 0 ',°’ P ' and K1: “mt-diumal tides (»=2I= 2, 3,4, 5): N2, M2, S2, 

equipotemitd w '1°^ ^ SUrfa “ ' opo * ra P 1 ’)’' “mplete to degree and order 15; (5) 

j i j | . * 09 or ^timeter biases; (6) correction to significant wave height H%n ( 7 } 

s r mu ,aser s,an r = <*> - <— 

set rarrc re on f re selecKd . saKiH,es> and 

satellites. Kaula s constraint equation [Kaula, 1966], which was inferred from surface gravitv 

’ were 5 ener atea. 1 hu- 1 b did not include direct altimeter data. 

5. ACCURACY EVALUATION 

“*! Cal ! bme *** accurac y of *e UT gravity models were performed, 
o panson of orbit fits using different gravity fields for Starlette, Ajisai, Seasat and Geosat were 
pe ormed. It is shown that using TEG-1, a Starlette five-day orbit fit is at the ~20 cm level, Ajisai 
ive- ay orbit fit is at the -15 cm level, and that a Seasat six-day orbit and a Geosat 17-day orbit have 


9 



Table 1. 

Satellite Data for the University of Texas 
Gravity Model, TEG-1 


Satellite 

Launch 

Data 

Inclination 

Eccentricity 

Altitude 


Date 




(km) 

Vanguard- 1 

1958 

Opticalt 

34° 

0.190 

2318 

Vanguard-2RB 

1959 

Opticalt 

33° 

0.183 

2318 

Courier- IB 

1965 

Opticalt 

to 

oo 

0 

0.016 

1100 

Geos-1 

1965 

Laser 

59° 

0.072 

1600 

BE-C 

1966 

Laser 

41° 

0.026 

1130 

DI-C 

1967 

Laser, Opticalt 

40° 

0.053 

1000 

DID 

1967 

Laser, Opticalt 

39° 

0.085 

1200 

Oscar- 14 

1967 

Doppler 

89° 

0.005 

1100 

Geos-2 

1968 

Laser 

o 

ON 

o 

0.033 

1400 

Peole 

1971 

Laser, Opticalt 

15° 

0.015 

650 

Geos-3 

1975 

Laser 

115° 

0.002 

830 

Starlette 

1975 

Laser 

0 

O 

0.020 

900 

Lageos 

1976 

Laser 

110° 

0.004 

5900 

Seasat 

1978 

Laser, Doppler, 

108° 

0.002 

800 



Altimeter and 






Crossover 




Nova-1 

1980 

Doppler 

vO 

O 

o 

0.002 

1200 

Geosat 

1985 

Doppler, Altimeter 






and Crossover 

108° 

0.000 

800 

Ajisai 

1986 

Laser 

50° 

0.001 

1500 

t Optical data currently withheld from gravity field solution 


Surface Gravity Data 

1° x 1° terrestrial mean gravity anomaly from 
Ohio State University [Rapp, 1986] 


crossover residuals at the ~25 cm level. Table 2 shows the summary for the Geosat orbit fits. Gravity 
field comparison using surface gravity data and a comparison of estimated TEG-1 ocean tidal 
parameters with solutions derived by other studies were also performed. Covariance matrices for 
TEG-1 and TEG- IS were calibrated to obtain estimates of errors associated with the gravity field 
using the consider covariance calibration technique [Yuan et al., 1988]. The predicted radial orbit 
errors using TEG-1 gravity field covariance matrix for Topex and Geosat are 13 cm and 24 cm, 
respectively (Table 3). 


6. CONCLUSION 

In this investigation, two gravity models, TEG-1 and TEG- IS, each complete to degree and order 
50 plus resonant coefficients, were generated. Ground-based tracking data collected by 14 satellites, 
altimeter crossover and surface gravity data were used to determine the TEG- IS gravity field model. 
TEG-1 contains Seasat and Geosat direct altimeter data in addition to all the data in TEG-1S. The 
gravity field models were derived simultaneously with orbit, ocean tides, quasi-stationary sea surface 
topography, and other geophysical parameters as well as the relative weights for each data set. The 
fields were evaluated using both data included and data withheld from the solution. Formal 
covariance matrices were calibrated to reflect realistic error estimates of the gravity field. Evaluations 
based on orbit fits and gravity anomaly residuals indicate that the gravity models have achieved a 
significant advancement over previously existing gravity models. 
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Table 2. Gravity Field Accuracy Evaluation ] 

Using Geosat Orbit Fits 


r oy v o» Cr* daily C D , density correction parameters adjusted 

Epoch 

TEG- IS 

TEG 1 

17-day orbits 

(rms) 

(rms) 

86/1 2/7 Doppler (cm/sec) 

0.67 

0.62 

Crossover (cm)t 

25 

22 

Altimeter (cm)t 

180 

32 

87/01/7 Doppler (cm/sec) 

0.62 

0.61 

Crossover (cin)t 

24 

25 

Altimeter (cm)t 

180 

32 


t Data types used for residual prediction only; altimeter data smoothed to 
represent gravity spectrum to (50x50) 


Table 3. Gravity Field Accuracy Evaluation 
Using Covariance Analysis 

Model 

Predicted Topex Radial 
Orbit Error (cm) 

Predicted Geosat Radial 
Orbit Error (cm) 

GEM-Tl 

25 

54 

TEG-1 

13 

24 

Topex orbit: 65° inclination, 1354 km altitude 
Geosat orbit: 108° inclination, 800 km altitude 
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